Numerical Modeling of Drying Residual RP-1 in Rocket Engines 

Alok Majumdar. Robert Polsgrove and Bruce Tiller 

Thermodynamics and Heat Transfer Group 
Marshall Space Flight Center. Huntsville, Alabama 


Introduction 

When a Rocket Engine shuts down under a fuel rich environment, a significant amount of 
unburned RP- 1 is trapped in the engine. It is necessary to clean the residual RP- 1 prior to 
subsequent firing to avoid any explosion due to detonation. The conventional method is 
to dry RP-1 with inert gas such as Nitrogen or Helium. It is difficult to estimate the 
drying time unless the engine is adequately equipped with instruments to measure the 
trace of RP-1 during the drying process. Such instrumentation in flight hardware is often 
impractical and costly. On the other hand numerical modeling of the drying process can 
provide a good insight for a satisfactory operation of the process. A numerical model can 
provide answer to questions such as a) how long it takes to dry, b) which fluid is a better 
dryer for RP-1, c) how to reduce drying time etc. 

The purpose of the present paper is to describe a numerical model of drying RP-1 trapped 
in a cavity with flowing nitrogen or helium. The numerical model assumes one- 
dimensional flow of drying fluid in contact with liquid pool of RP-1 (Figure -1). An 
evaporative mass transfer takes place across the contact surface. 



Figure -1 Drying of liquid RP- 1 by gaseous Nitrogen 

The considered problem was to find the rate of mass transfer of RP-1 into nitrogen stream 
for a given pressure, temperature and flowrate assuming no heat transfer between RP-1 
and nitrogen. It is also assumed that mass transfer between liquid and gas occurs through 
a thin film separating two streams. This film contains saturated vapor at interface 
temperature, which is assumed to be at liquid temperature. 



Governing Equations 


The equations governing the mass transfer of liquid to a gaseous stream can be expressed 
as 1 : 

Mass Transfer Rate. Wa (lb- mo 1/sec) 
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Where x A o is the mole fraction of vapor of the saturated film near the liquid surface and 
x A ~ is the mole fraction of vapor in the gaseous stream. A is the surface area and k xm is 
mass transfer coefficient given by: 

Mass Transfer Coefficient, k m (lb- mol/ft 2 -sec) 




P mot Dab 


2.0+0.60 


P/V^L 

Pi 


. 0 . 5 , 


Pi 

PfD/ 


\ 0.33 


(2) 

L is a characteristic length scale, Pf and Lif are density and viscosity of vapor at film 
temperature. p mo | is the molar density, v« is the free stream velocity of gas near the liquid 
surface. The diffusivity between liquid and gas, Dab is given by: 

Mass Diffusivitv. Par (fr/sec) 
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The subscript A and B refer to liquid and gas respectively. The subscript c refers to 
critical point and p and T are pressure and temperature respectively. 
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RP-1 


M = 0.0745 lbm; V = 2.5 in 3 ; A = 5 in 2 ; T = 60 ° F 


Figure 2 - GFSSP Model of mass transfer 

GFSSP Model 

A simple transient GFSSP (Generalized Fluid System Simulation Program") model was 
constructed to estimate the evaporative mass transfer rate ot RP-1 using nitrogen and 
helium as drying fluid. The model is shown in Figure 2. Figure shows that nitrogen is 
flowing through a duct and is in contact with RP-1 stored in a cavity. Node 1 and 3 are 
boundary nodes where pressure and temperature are specified. Node 2 and 5 are internal 
nodes where all scalar properties such as pressure, temperature and concentrations are 
calculated. Node 5 contains RP- 1 of known initial mass, volume and temperature. 

The mass transfer model described in equations 1 to 3 was implemented in User 
Subroutine of GFSSP. User Subroutine is an effective way of incorporating new physical 
model in the code. The GFSSP process flow diagram is shown in Figure 3. The code 
consists of tliree major modules: 1) Preprocessor, 2) Solver and Property Program and 3) 
User Subroutines. User creates the input data file with the help ot the Preprocessor. 
GFSSP has two preprocessors: command line and visual. Solver and property program 
generates and solves all pertinent equations with the help of fluid property routines. User 
subroutines are set of blank subroutines called from solver module. This allows the user 
an opportunity to add or modify the governing equations according to their need. The 
developed user subroutine is then compiled and linked with solver module to create a 
customized executable. 
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Figure 3. GFSSP Process Flow Diagram 


Results 

Several cases were run to study the effect of various parameters on drying time. The 
results of the parametric were shown in Table 1. Results indicate that a) Nitrogen is a 
better dryer of RP- 1 than helium, b) Drying rate significantly increases at higher RP-1 
temperature and c) Lower pressure and higher flowrate increases drying process. 
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Table 1 . Results of Parametric Study 
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8.59 

1.27 
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7.93 

0.85 
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13.59 

1.27 
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13.59 

1.26 

55 
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100 
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0.183 

5.86 

1.63 

25 

20 
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60 

0.039 

8.6 

0.196 

25 

1 

100 

60 

0.833 

12.4 

0.238 


3.146 

5.001 

1.922 

I. 935 
3.523 
20.002 

II. 67 

























































